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SUMMARY 


The  problem  of  the  prediction  of  the  exact  pressure  field  produced 
at  close  ranges  by  an  impulsive  sound  source  near  a  fluid-layered  solid 
interface  is  considered.  A  theory  is  developed  that  should  accurately 
provide  arrival  time  data  and  exact  pressure  response  waveform  data  for 
a  model  consisting  of  an  explosive  acoustic  source  in  a  fluid  overlying 
a  two-layered  solid.  Experimental  pressure  response  curves  have  been 
obtained  with  the  use  of  models  and  are  presented  to  support  the  theo¬ 
retical  development.  Although  the  final  results  of  the  theory  are  not 
at  present  in  closed  form,  it  is  anticipated  that  such  results  should 
be  possible  in  the  near  future. 

The  method  of  Cagniard,  which  is  a  Laplace  transform  technique,  is 
utilized  in  the  mathematical  development.  Although  this  type  of  deri¬ 
vation  becomes  quite  complicated,  it  is  worthwhile  because  it  enables 
one  to  obtain  expressions  describing  exact  pressure  response  waveforms 
that  are  in  closed  form;  that  is  requiring  no  further  integrations. 

It  was  shown  in  a  previous  paper^  that  for  low-frequency  sound 
transmission  near  the  ocean  floor,  the  bottom  cannot  be  considered  as 
a  liquid  because  of  the  possibility  of  surface  waves  occurring  at  the 
bottom  interface.  These  waves  are  attenuated  in  a  direction  normal  to 
the  interface,  and  hence  are  not  important  at  high  frequencies  or  large 
distances  from  the  interface. 

The  experimental  data  presented  was  obtained  in  a  model  using 
previously  developed  experimental  techniques.^  The  sound  source  was  a 
coaxial  spark  gap  and  the  hydrophone,  or  detector,  was  a  small  ultrasonic 
probe.  Model  solids  included  plexiglas,  plaster,  flint  glass,  aluminum, 
and  lead. 

It  is  concluded  that  the  theoretical  approach  of  Cagniard  appears 
to  be  quite  applicable  to  problems  of  the  type  presented  here.  The 
extended  theory  developed  in  this  report  is  accurate  in  predicting  most 
expected  arrivals  at  short  ranges  in  the  experimental  model.  It  is 
anticipated  that  further  theoretical  development  will  provide  exact 
waveform  data.  Except  for  rare  situations,  it  is  found  that  surface 
waves  cannot  exist  at  real  solid-solid  interfaces.  Additional  surface 
waves  appear  to  exist  at  the  fluid-solid  interface  because  of  reflections 
from  lower  solid-solid  interfaces.  The  results  obtained  in  this 
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investigation  continue  to  illustrate  the  large  errors  in  wave  propagation 
predictions  near  an  interface  that  can  occur  when  the  ocean  bottom  is 
considered  to  be  a  fluid.  It  is  hoped  that  further  work  on  problems  of 
this  type  will  eventually  permit  the  prediction  of  the  exact  pressure 
field  produced  by  an  explosive  source  near  a  real  ocean  floor. 
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INTRODUCTION 


A  theory  recently  published  by  Strick^  has  been  found  to  give 
reasonable  accuracy  in  prediction  of  the  pressure  fields  at  short  ranges 
from  an  impulsive  sound  source  near  a  fluid-homogeneous  solid  interface. 
However,  the  assumption  of  the  homogeneous  solid  greatly  restricts  the 
application  of  this  theory.  When  considering  low  frequencies  it  is 
difficult,  if  not  impossible  to  find  an  ocean  bottom  that  is  homogeneous 
down  to  a  depth  of  at  least  one  wavelength,  the  point  at  which  surface 
wave  effects  can  usually  be  neglected.  If  the  mean  frequency  under  con¬ 
sideration  is  15  cycles,  a  typical  hard  bottom  (a  =  7500  ft/sec)  would 
have  to  be  homogeneous  down  to  a  depth  of  about  500  feet. 

Hence,  this  report  is  intended  to  extend  the  theory  advanced  by 
Strick  to  a  model  consisting  of  a  fluid  overlying  a  two- layer  solid. 
Although  the  final  result  of  the  extended  theory  is  not  in  closed  form, 
it  is  developed  to  a  point  that  one  can  predict  and  identify  the  elastic 
wave  arrivals  occurring  in  the  pressure  response  curves  obtained  from 
the  experimental  model. 

Pekeris^  presented  the  normal  mode  theory  for  explosive  sound  trans¬ 
mission  in  three  liquid  layers  in  1948,  and  this  problem  has  been  further 
investigated  by  Officer‘d  and  Press  and  Ewing^.  McLeroy®  has  published  a 
complex  image  theory  for  near- fie Id,  low  frequency  propagation  in  shallow 
water  overlying  a  many- layered  bottom,  but  it  is  not  valid  for  the  problem 
to  be  considered  here. 


MATHEMATICAL  PRESSURE  RESPONSE  DERIVATION* 
FOR  THE  TWO-LAYER  BOTTOM  MODEL 


The  geometry  of  the  problem  is  presented  in  Figure  1.  As  in 
Reference  1,  the  liquid  is  considered  to  be  semi-infinite  in  extent. 


*To  the  best  of  our  knowledge  this  derivation  has  not  been  made  elsewhere. 
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FIGURE  1.  THE  ASSUMED  MODEL 


It  will  be  seen  that  the  derivation  for  this  extended  problem  will 
closely  follow  that  given  by  Strick^  in  his  original  derivation  for  a 
liquid-homogeneous  solid  model.  A  simplified  version  of  this  derivation, 
along  with  a  number  of  added  explanatory  comments,  can  be  found  in 
Reference  (1).  It  will  be  assumed  that  the  reader  is  familiar  with  at 
least  one  of  the  aforementioned  papers. 

Figure  2  serves  to  identify  the  11  displacement  potentials  that 
will  be  considered  in  this  problem.  It  will  be  assumed  that  any  waves 
existing  after  two  reflections  from  the  lower  (solid-solid)  interface 
can  be  neglected.  For  many  situations  this  is  not  a  good  assumption, 
but  it  is  not  so  unreasonable  when  one  realizes  that  the  predominant 
waveform,  the  surface  wave,  dies  out  exponentially  with  distance  normal 
to  the  Interface.  The  introduction  of  the  middle  layer  into  the  problem 
has  made  a  rigorous  solution  of  the  problem  by  the  potential  technique 
difficult,  as  there  are  now  an  infinite  number  of  potentials  that  should 
be  considered.  As  in  the  homogeneous  bottom  problem,  the  pressure 
response  will  be  obtained  for  a  delta-excited  line  source,  instead  of  an 
explosive  point  source. 

cp^  is  the  displacement  potential  due  to  the  line  sound  source  in  an 

infinite  fluid,  cpf  represents  the  displacement  potential  due  to  energy 
propagated  back  inro  the  liquid  because  of  the  presence  of  the  fluid- 

solid  interface.  Likewise,  cpf  is  the  displacement  potential  representing 
energy  propagated  back  into  the  fluid  because  of  the  presence  of  the 
solid-solid  interface.  As  shown  in  Appendix  A  of  Reference  1,  two  ' 
potentials  are  necessary  to  describe  the  wave  propagation  in  an  infinite 
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FIGURE  2,  THE  11  DISPLACEMENT  POTENTIALS 

Si 

solid  tnediuni.  Hence,  potentials  such  as  cp^  and  are  necessary  in 
the  solid. 

Tbs  displacement  potential  in  an  infinite  fluid  due  to  a  harmonic 
line  source  is  given  by  ^ 


f,  ia)t. 
cp^(x,z;e  ) 


^-iuj(rz  +  qx) 


where 


P 

o 


Zip^o)^ 


and 


q 


-  sin  0 
c 


(1) 


-II- 


the  angle  0  is  defined  as  the  angle  between  the  vertical  z-axis 
and  the  wave  vector  associated  with  each  elementary  plane  wavefront 
assumed  in  obtaining  Equation  (1).  The  other  potentials  in  the  problem 
are  given  by  Equations  (2a-j) ,  where  A(q)  through  N(q)  represent  re¬ 
flection  and  transmission  coefficients  which  will  be  evaluated  shortly:. 


(2a) 


doe 


tu)t 


(2b) 


\|i^  (x,z;e  )  =  Qj 


1  2  ^ 

-iu)[(^  -  q  )  z  +  qx]^ujt 


C(q)e  " 


(2c) 


1.  iUDtv 

cp«  (x,z;e  ) 


=  Q. 


-icu[(-y  -  q^)  (2H  +  2J  -  z)  +qx] 

D(Q)e  a  _ dae  (2d) 


1  n  ^ 

-iU)[(-*T  -  q  )  (2H  +  2J  -  z)  +qx] 


=  Q, 


E(q)e  b 


-iae 


iu)t 


(2e) 


CP3  (x,z;e  )  = 


-i(i)[(^  -  q^)  (2J  +  z)  +  qx] 

F(q)e  a^ _ dqe^^*^ 


(2f) 


,  1,  itut.  - 

(x,z;e  )  = 


■1U3[(-^  -  q^)  (2J  +  z)  +  qx] 


G(q)e  b 


dqe 


i(J3t 


(2g 


2.  iuut. 
cpj^  (x,z;e  )  = 


1  2 

“  -iu)[(-2  -  q  )  z  +  qx] 

r  M(q)e  _ da^ 


(2h) 
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(21) 


(x,z;e  )  =  Q 


1  2  ^ 

■iutr  (—7  -  q  )  z  +  qx] 


I  Cq)  e  k 


dqe 


i'i)t 


cp2(x,z;e  )  = 


IJCq)e 


-Lii)[r(2H  +  2J  -  z)  +  qx] 


dqg. 


1(1)  t 


r 


(2j) 


The  boundary  conditions  are: 

(1)  The  normal  stresses  are  continuous  at  both  interfaces. 

(2)  The  normal  displacements  are  continuous  at  both  interfaces. 

(3)  The  tangential  stresses  are  zero  at  the  liquid-solid  interface 
and  continuous  at  the  solid-solid  interface. 

(4)  The  tangential  displacements  are  continuous  at  the  solid-solid 
interface . 

These  boundary  conditions  give  10  boundary  equations  which  can  be 
used  to  evaluate  the  10  unknowns,  A,  B,  C,  D,  E,  F,  G,  M,  I,  and  N. 
Actually,  as  we  are  only  interested  in  the  response  in  the  fluid,  it  is 
only  necessary  to  determine  A(q)  and  N(q) ,  The  details  of  the  derivation 
of  these  quantities  are  given  in  Appendix  A.  As  expected,  A(q) ,  or  A(u) 
with  a  change  of  variables,  is  unchanged  from  the  similar  reflection 
coefficient  occurring  in  the  original  fluid  -  homogeneous  solid  derivation 
and  is  given  by 


A(u)  = 


Y(K2  -  2u^)  +  4u^Q'yP  -  ^  K2  a 

2 

Y(Kt  -  2u^)  +  Au^c^yP  +  *^7  O' 

t  P  S  ^  4 


where 

2  2  ^ 

a  =  (Kj^  -  u  )  , 

9  2  ^ 

P  =  (K2  -  u  )  , 


(3) 
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and 


Y  =  (1  -  u  )  , 

Kj  =  c/a  , 


K2  =  c/b  , 


u  =  cq 

c  =  compressional  sound  velocity  in  the  fluid, 
a  =  compressional  sound  velocity  in  the  first  solid, 
b  =  shear  sound  velocity  in  the  first  solid. 


N(q)  is  considerably  more  complicated  than  A(q)  and  is  given  by 

-2i.J(a'  -  D  -2ia)J(b'  -  O  -ia.J(a'  +  b'  -  2r) 

N(q)  - 3 - ^ 

3®  4^ 


where:  ^  2 

w.pf +  4p“ibVi.'r+  P  't-d 


(5) 


=  [(1  -  2b^q^)  (4p^r)(H^e^  -  H^e^)]  -5-  [a'  -^  +  Tn] 


£f_ 


N2  =  [16  p^rb^^a'b'  (^223  -  ^2^3^^  ^ 


(6) 


(7) 


N3  =  [-8pWq(l  -  2b^q^)]  [b'(H^e3  -  e  ^H3)  +  ^ 


[a'  ^  +  rn] 


ps 


(8) 


C3  =  {[2p  Vqa'  -  ^  (1  -  2k^q^)]  +  [2k^q^  +  (1  .  2b^q^)  ]  • 

®2  ^  ^ 

[-2p  Vg'q  +  ^  (1  -  2k^q^]}  (9 
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2.  » 


{[p  ^(1  -  2b V)  -  (1  -  2k^q^)]  +  [-2kV'  + 


[2p  \^g'q  -  (1  -  2k^q^)]| 
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C2p  ^k^g'q  -  (1  -  2k^q^)]| 


{(^  -  k^ 

{fq 


2k^q^  +  ^  (1  -  2b^q^ 

®2 

P 


2k^qb'  -  b^b'q]]- 


>1 


{(•'  -  k^<«'  -^  ^)[^;  «  -  + "''“]} 


{(’  ■^  M)  G' 


p 

(10) 

]  +  2k^^] 

(11) 

+  2k^qb ■ 

(12) 

(13) 

(14) 

(15) 

(16) 
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(17) 


g  =  cotnpressional  sound  velocity  in  the  lower  solid 

k  =  shear  sound  velocity  in  the  lower  solid 

2 

n  =  (1  -  2b^q^)  +  4b^q^a'b'  (18) 

The  total  displacement  potential  in  the  fluid  is  given  by 


f  f  ^  f  ^  f 

9  +  9i  +  ^2  . 

and  the  total  fluid  pressure  response  is 

f  f  ■  f  .  f 

P  =  Po  +  Pi  +  P2  • 

cp^  ,  and  hence  p^  ,  is  unchanged  from  the  original  homogeneous 
solid  problem.  For  an  impulsive  type  sound  source. 


,  P  c  H(t  -  — ) 

f .  ...  o  c 

p^(x,z;  delta)  =  - 2~2 -  h 


(19) 


0^2  V 


where 

2  2  ^ 

r^  =  (x'"  +  z^)  , 

H  =  the  unit  step- function, 

and  P  is  an  undetermined  constant  depending  on  the  strength  of  the  sound 
o  f  f 

source.  Likewise,  since  cp^  is  unchanged,  p^  is  unchanged,  and  is  given  by 
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where 


f(u)  =  (kI  -  2u^)  -  4u^(u^  -  K^)  (u^ 

r,  r 

(c)  Case  1:  t  s  -i  Case  2:  t  — 
c  c 


p^(x,z;  delta) 
where 

D(u)  =  yfiu)  +  (^)  a  , 

It  is  now  necessary  to  find  p2  (x,zj  delta),  the  pressure  response 

in  the  fluid  due  to  contributions  from  the  soUd-sol  d  boundary  when 
the  sound  source  is  impulsive.  It  was  shown  above  that  the  integral 

expression  for  p^  is 


P  c 
o 


Kr  |f(u)i  -  k“  |oi 

P 


(21) 


f 

P2 


(x,z;e  ) 


Pq  ’’  +  2J  -  z)  +  qx]  j^^luut 

n  J  f 

-CO 


(22) 


Or,  using  Equation  (4) 


r.  lujt. 

P2(x,z;e  )  == 


2i 


N  gi^[r(2H  -  z)  +  2Ja'  + 


rw<=^.3  -  53.3^) 


®  -ioD[r(2H  -  z)  +  2Jb'  +  qx]  iout 


2i 


2i 


P  N^e 


Mqe" 


rW(H  C  -He) 
4  3  3  4 


fj  g*iuj[r(2H  -  z)  +  J(a'  +  b')  +  qx]^^^iu)t 


r«(V3  -  -3*4) 
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=  [p2(x,z; 


-I  ,  r  f ,  iUUtT 


+  [P2(’«.z»  (23) 


The  expression  for  the  s-multiplied  Laplace  transform  of  the  p^ 

contribution  to  the  total  pressure  response  due  to  a  step-function 
excitation  of  the  source  can  now  be  obtained  by  replacing  iti)  by  s  and 

dropping  the  e  in  Equation  (23).  This  is  a  technique  used  in  filter 
analysis  and  brief  explanations  can  be  found  e Isewhere , ^ 


Thus , 


f  ,  i^^t  _  o 
P2  )  -  2i 


®  N  p"®Cr(2H  -  z)  +  2ja'  +  qx] 
P  P  1 


dq 


rW(H^e3  -  H3 


"if 


»  „  -s[r(2H  -  z)  +  2Jb'  +  qx] 

p  N,e  ^ 


dq 


rw(3j.3  -  =3,^) 


-if 


N  -  z)  +  j(a'  +  b')  +  qx]^^ 

■  ■=■3*4) 


(24) 


The  s-multiplied  Laplace  transform  of  a  sum  is  the  sum  of  the  transforms 
of  the  individual  components,  so  it  is  possible  to  separate  the  problem 
into: 

P2(x,z;step)  =[p2(x,z;step) ]j^  +  [p2(x,z;step)]2  +  [p2(x,z ; step) , 


where 


-  -[(2H  -  z)  Y  +  2c Ja'  +  ux] 

®  N-e  du 

fx  _ o 


21  J 


r“<S4-3  -  H3  .,) 


-  ~[(2H  -  z)  y  2cJb'  +  ux] 
—  P  “  N,e  ^  du 

(P2>2  =  2l 


rw(E^.3  -  =3.,) 


-19- 


CO 

p 


(25) 


/  fv  O 

^^2^3  2i 


N^e 


"  [(2H  -  z)  y  -h  cJ(a*  +  b’)  +  ux] 


du 


™<“4*3  -  "S‘4) 


and 


U  =  cq,  r  =  . 

To  obtain  the  inverse  s-multiplied  Laplace  transform  it  will  be 
necessary  to  have  an  expression  of  the  form 


P(x,z; 

Then,  using 


step)  =  (Constant) 


the  derivative  rule 


00 

[fn(u)]e  ^^dT  , 

0 

for  s-multiplied  Laplace 


(26) 

transforms, 


j:(f*(t))  =  s  T(f(t))  -  s  f(0)  , 

and  the  fact  that  the  delta-response  is  the  derivative  of  the  step-response, 
it  will  be  possible  to  obtain  the  delta-function  excitation  pressure  response 
in  a  closed  form: 


p(x,z;  delta)  =  Constant  [fn(u)] 


So,  to  reduce  Equation  (25)  to  three  equations  of  the  form  of  Equation  (26) 
it  is  necessary  to  define  t^,  and 


T 


T 


1 


2 


3 


(2H  -  z) 
c 

(2H  -  z) 
c 

(2H  -  z) 
c 


9  %  1  2  % 

u  )  +  2J(—  -  — )  +  — 

a  c 


o  %  .  2  % 

u  )  +  2J(-  -  — )  +  - 

b  c 


DC  a  c 


2  % 


1 


2  h 


(27a) 

(27b) 

(27c) 
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'X T 2  1  each  contain  at  least  one  term  that  is  not  present 

in  the  transformation  to  t  in  the  homogeneous  bottom  problem.  The 
transformation  in  that  case  is 

h 

T  (Homogeneous  bottom)  =  ~  -g)  q  _ 

X  c  c 

It  was  relatively  easy  to  determine  u  as  a  function  of  T  in  that  case, 

but  for  the  current  problem  this  becomes  much  more  difficult.  On  the 

other  hand,  it  is  still  possible  to  obtain  an  expression  for  du  as  a 

function  of  dr.  and  u. 

1 


From  Equations  (27a),  (27b)  and  (27c),  after  a  substitution  of  previously 
defined  quantities,  it  is  found  that: 


dT^ 


du  = 


du  = 


du  = 


X 


a  c^  c 


-(2H  -  2)u  2Ju  .  X 

— * - -  -  7  +  — 

yc  '' 


dT, 


0 


‘(2H  -  z)u  2Ju  X 

“  ^  c 


yc 


- 


dr. 


""(2H  -  z)u  Ju 
yc  2 


(28a) 


(28b) 


(28c) 


Then, 


P2(x,z;  step)  =  -~ 


-ST 

N^e 


yWv-j^  ^“4®  3  "  “3®4^ 


-st2 

N2e  dT2 


YWV2  (2^, 


4*3  *^3  '4 


«/.)  2i 


-ST 

^ 


YWV3  (-4*3  -  "3*4) 


(29) 


Tz  '3 

where  is  defined  by  the  denominators  of  Equations  (28a),  (28b)  and  (28c). 
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The  problem  has  been  reduced  to  expressions  of  the  form 


step)J  =  (Constant)  T  fiu(u)Je 


T. 

1 


(i  =  1,  2,  3) 

But,  since 

s[p2(x,z;  step)]  =  P^(x,z;  delta) 
as  was  discussed  previously. 


j^p^CxjZ;  delta)J  =  s(Constant)  ^^^i 

^  c 

T  . 

L 

(i  =  1,  2,  3)  (30) 

The  right-hand  side  of  this  expression  is  in  the  form  of  an 
s-multiplied  Laplace  transform,  except  for  the  fact  that  the  limits  on 
the  integral  are  not  proper.  Once  the  integral  over  the  contour  c 

i 

is  changed  to  an  integral  over  =  0  to  expression  (30)  will  be 

of  suitable  form  to  take  the  inverse  Laplace  transform.  The  problem  will 
then  be  in  a  closed  form  requiring  no  further  integration.  This  transform 
to  a  proper  form  of  integral  is  complicated  by  the  presence  of  a  large 
number  of  poles  and  branch  points,  and  by  the  fact  that  u  cannot  be  found 
as  an  explicit  function  of  t,.  Work  is  continuing  at  the  present  time  on 
this  transformation.  But,  tftere  is  still  a  large  amount  of  extremely 
useful  information  in  the  present  integral  form  of  the  delta  response. 

The  closed  form  solution  is  desired  to  theoretically  predict  the  exact 
waveform  of  the  various  elastic  waves,  but  this  waveform  can  also  be  found 
experimentally  in  a  model.  Thi  experimental  problem  then  reduces  to 
identifying  the  type  of  wave  with  the  waveform.  The  present  form  of  the 
theory  is  sufficient  to  make  this  identification. 
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DISCUSSION  OF  THEORETICAL  RESULTS 


In  the  fluid  -  homogeneous  bottom  problem  Strick  found  that  each 
elastic  wave  could  be  predicted  by  the  existence  of  poles  and  branch 
points  in  his  expression  for  the  total  pressure  response.  Thus,  the 
indicated  procedure  at  this  point  is  to  search  for  similar  singularities 

in  Equation  (29) ,  noting  that  in  general  p^  and  p^  have  already  pre¬ 
dicted  the  existence  of  P,  S,  direct,  reflected,  Rayleigh  and  Stoneley 
waves  due  to  the  presence  of  the  water-solid  interface . ^ ^ 


Branch  points  are  present  in  the  denominator  of  Equation  (29)  at 
c  c  c  c 

^  “  -  ^9  “  ~  k**  When  these  values  for  u  are  inserted  into 

the  expressions  for  t^,  it  is  found  that: 


(31a) 


r,(u  =  f) 


Ti(u  -  1) 


^  L2H  -  z) 
c 

^  C2.H  -  2) 

c  ^ 

=  (2H  -  z) 

c  ^ 


+2J(^ 

b  a 


+  2J(^ 


1 

—y)  +  2  J 


-P1S2P2R  P^s^P 


p  p  p  p  p 

r2^3  2^1 


1 

(31b) 

(31c) 

(3  Id) 


In  the  transformation  from  u  to  the  y  factor  should  drop  out  and 


be 


C^tV 

replaced  by  6  1  -  — »  representing  a  reflected  wave  from 


1 


the  liquid- solid  interface.  See  the  original  paper  by  Strick  for  an 
elaboration  on  this  point. 
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T2(u  =  ;)  =  (1  -  ^)  +  2J(^  -  -1)  +  I  t  3 

a  ba  1  Z  z  Z  Z  1 


/  c.  (2H  -  z)  c  V  .  X  - 

~  c  b  P^S^Pj^ 


Ole) 

(31f) 


^2^^^  -  g>  ■  ^  ^  c  ■  2^  ^'^^  2  '  2^  g  ■■  "P  S  P  S  P 

®  g  h  g  1  2  J  Z  1 


T  (u  =  h  =  (1  -  ■^)  +  2J(-j  -  -j)  +f^t  (31h) 

2  k  c  1^2  k  P^b2b3b2P^ 


1  1.  ,  X 


T3(u  =f)  = 


i- 

iln  -Z).  Cv  .T/l  In  .X 

^  (1  -  -)  ^  J(  -  2)  •"  a 

a  b  a 


(311) 


c.  (2H  -  z 

,3(u  =  -)  =  ^ 


(1  -  — )  H-  J(—  -  — )  +  — 
b  a  b 


(31j) 


(„  -  f)  -  (1  -  %)  +  j[(4  -  ^)  +  (^  -  ^)  ]  +  f 

®  g  b  g  a  g  ° 


“  '•'l"2V2''l  ■  ‘‘'l‘'2‘'3=2‘’l 


’3<“  -  k>  ■  -  ,2>  ^  1<^2  .2)  ,2)  J  +  k 

“  'PiP2S3"2‘'i  ■  '•'lS2V2''l 


(3U) 
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Appendix  B  contains  geometrical  derivations  for  the  travel  times 
of  the  critically  refracted  and  the  reflected  -  critically  refracted 
arrivals  that  are  to  be  expected  in  this  problem.  An  analysis  of  these 
travel  time  expressions  will  show  that  the  geometrical  travel  times  are 
given  by  equations  identical  to  those  given  by  Equations  (31),  except 
that  no  geometrical  arrivals  can  be  found  to  correspond  to 
Equations  (31i)  and  (31j).  Also,  an  analysis  of  the  experimental  re¬ 
sults  will  show  that  no  arrivals  are  present  corresponding  to  those 
described  by  (31i)  or  (31j).  Note  that  when  (31i)  is  not  complex, 

(31j)  is  complex,  and  when  (31i)  is  complex,  (31j)  is  not  complex. 
Hence,  arrivals  corresponding  to  (31i)  and  (31j)  could  never  exist 
simultaneously. 


Poles  exist  in  Equation  (29)  corresponding  to  the  vanishing  of 
W,  or  (^^€^  “  ‘  ^1^  ^2  given  by  the 

denominators  of  Equations  (28a,  b,  c)  respectively.  It  is  hypothesized 
that  when  the  value  of  u  that  corresponds  to  the  vanishing  of  is 

substituted  into  the  corresponding  equation,  the  equation  will 

give  the  travel  time  of  an  arrival  that  has  been  refracted  at  the  fluid- 
solid  interface,  reflected  at  the  solid-solid  interface,  and  again 
refracted  at  the  fluid-solid  interface.  Geometrically,  the  approach 
necessary  to  determine  a  travel  time  for  such  a  reflected  wave  consists 
of  using  an  iterative  process  to  find  a  reflected  travel  path  that 
satisfies  Snell's  law  for  a  given  set  of  dimensional  and  elastic  parame¬ 
ters.  It  will  be  seen  that  the  two  methods  give  identical  results,  so 
it  is  concluded  that  the  mathematical  approach  hypothesized  is  valid. 

From  Equation  (5) , 


W  =  p  +4  b  q  b  T  +  p  (1  -  2b  q  )  — 

P  ^ 

=  [y(K2  -  2u^)  +  o'] 


(32) 


The  term  in  brackets  in  Equation  (32)  is  exactly  the  denominator  in 
Equation  (3).  In  the  liquid-homogeneous  solid  problem  the  vanishing  of 
the  denominator  of  Equation  (3)  defined  the  phase  velocity  of  the  Stoneley 
interface  wave.  Hence,  Equation  (32)  simply  predicts  a  Stoneley  wave 
traveling  along  the  fluid-solid  boundary.  Using  the  data  from  Table  1 
and  the  Stoneley  wave  charts',  it  is  found  that  for  a  water-plexiglas  inter¬ 
face,  Equation  (32)  vanishes  for  a  value  of  u  corresponding  to  a  Stoneley 
velocity  cf  1.03  x  10^  centimeters  per  second. 
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The  expression  -  H^e^)  contains  terms  depending  on  the 

elastic  constants  of  both  solid  media,  and  it  is  hypothesized  that  the 
vanishing  of  this  term  describes  the  propagation  of  a  Stoneley  surface 
wave  at  the  interface  of  the  two  solid  media.  It  seems  likely  that  the 
expression  -  H^e^)  should  reduce  to  the  equation  obtained  by 

Stoneley  to  describe  the  propagation  of  a  surface  wave  at  the  inter- 
face  of  two  solid  media , 

If  multiplied  by  -  kq  p  and  the  resulting 

expression  is  simplified,  it  is  found  that 

-k'c'S^^(E^e^  -  = 


+  (p'l)"  c'^g'k'  -f 


S,  2 


s«  2  s. 


+  (p‘^)‘  -  2p”V^  +  <p“l,'] 
+  c'^[4p  ‘p  ^b^c’^a’b'  -  4(p  k^c'^a'b'  -  4(p**)  b^c'^g'k' 


s.  s. 


S.  2 


+  4p'lp V^'2g.k'  +  4p'lpV  +  4p'^p  V  -  4(p'l)  V  -  4(p'2) Vj 


+  [c'^a'b'g'k’r4(p  h  -  8p  ^p  +  4(p^^)  +  4(p^b^’ Va’b 


-8p  ^p  ^b^c'Va'b'  +  4(p  k^c'^a'b'  +  4(p  k^c'^g'k' 

®1  ^  A  0  9  0  O  S«  2  y 

+  4(p  b  c'  g'k'  -  8p  ^p  ^  k^b^c'^g'k'  +  4(p  k^ 


2  s.  2 

,  +  4 
y 

s«  2 


s.  2 


4(p“b  V  -  8p'lp^Vb2], 


where  c*  =  “  . 

q 


(33) 


Now,  let 

a’  =  ±  ^  A^,  b'  =  ±  ^  B^,  g'  =  ±  A2,  and  k'  =  ±  . 
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Then 


k'c' 


'^p  ^*'"4*3  ~  "3*4^  "  (P  ■  P  "  (p  ^^2  ’’’  P  P  i 


s„  2  s,, 


2„  .1 


2r®l  ®9  ®on  o 

+  2Kc'  [^p  -  p  -  p  ^  +  p 

^12  ^2 

where  K  =  2(p  b  -  p  k')  . 


-  DCA^B^  -  1)J 


(34) 


This  is  exactly  the  equation  describing  the  Stoneley  wave  that  was 
obtained  by  Stoneley  in  his  original  solid-solid  problem.^ 

Scholte®  has  investigated  Equation  (34)  in  an  attempt  to  determine 
the  conditions  under  which  Stoneley  waves  can  exist  at  a  solid-solid 
interface.  He  found  that  they  can  exist  at  the  interface  only  under  a 
set  of  extremely  stringent  conditions.  The  results  of  Scholte's  investi¬ 
gation  are  presented  in  Figures  3  and  4.  Stoneley  waves  can  exist  only 
in  the  shaded  regions,  and  are  the  Lame ^  constants  of  the  solids. 

These  conditions  are  in  fact  so  stringent  that  it  is  difficult  to  con¬ 
ceive  of  a  model  that  can  satisfy  them.  It  is  necessary  to  consider  a 
radical  model  like  plexiglas  over  lead  in  order  to  formulate  a  system 
that  will  contain  Stoneley  arrivals  that  originate  at  a  solid-solid 
interface.  When  such  a  model  is  considered  experimentally,  the  Stoneley 
velocity  is  so  slow  that  multiple  reflections  occurring  at  the  boundaries 
of  the  experimental  system  makes  the  identification  of  a  ”Stoneley-type*' 
wave  impossible.  Certainly,  a  model  such  as  plexiglas  over  lead,  where 
the  densities  differ  by  a  factor  of  ten,  should  not  be  found  in  nature. 


EXPERIMENTAL  RESULTS 


The  apparatus  used  in  the  experimental  phase  of  this  project  was 
quite  similar  to  that  utilized  in  previous  work  with  the  homogeneous 
bottom.  A  brief  description  of  the  instrumentation  and  model  facility 
can  be  found  in  Appendix  C  of  this  paper.  The  major  experimental  problem 
encountered  in  the  layered  bottom  extension  was  a  question  concerning 
the  bonding  of  the  solid  layers.  As  described  in  Appendix  C,  it  was 
found  that  for  our  very  smooth  and  flat  surfaces  the  bonding  of  the 
layers  was  not  important.  The  two  models  used  in  the  experimental  phase 
of  the  project  were  water  -  plexiglas  -  plaster  and  water  -  glass  - 
aluminum.  The  elastic  and  physical  parameters  of  these  materials  are 
given  in  Table  1. 


(Text  Continued  on  Page  29) 
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3 


Asymptote  of  B 


FIGURE  3. 


CONDITIONS  UNDER  WHICH  STONELEY  WAVES  CAN  EXIST  AT  A 
SOLID- SOLID  INTERFACE  WITH  ^  ^  =  1 
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TABLE  1 


ELASTIC  PARAMETERS  IN  THE  MODEL 


Material 

Compressional  Velocity 
(cm/sec) 

Shear  Velocity 
(cm/ sec) 

Density 
(gm/  cm^) 

Thickness 

(cm) 

Water 

1,46  X  10^ 

- 

1.00 

“  (20.0) 

Plexiglas 

2.62  X  10^ 

1.30  X 

10^ 

1.  17 

3.81 

Plaster 

2.  83  X  10^ 

1.80  X 

10^ 

1.89 

®  (24.0) 

Glass 

3.30  X  10^ 

1.56  X 

10^ 

3.50 

1.27 

Aluminum 

5.26  X  10^ 

2.95  X 

10^ 

2.  70 

“  (20.9) 

Figures  5  and  6  are  typical  pressure  response  versus  time  curves 
obtained  for  the  two  models.  The  amplitudes  are  relative  and  travel 
times  are  accurate  to  ±  1.0  microsecond.  The  noted  wave  identifications 
will  be  discussed  in  the  next  section. 


EXPERIMENTAL  -  THEORETICAL  COMPARISON  AND  DISCUSSION 


The  large  number  of  elastic  wave  arrivals  occurring  in  a  short  time 
period  in  this  problem  prohibit  all  of  the  theoretical  arrivals  to  be 
identified  from  any  one  experimental  photograph.  In  an  attempt  to  ex¬ 
perimentally  identify  as  many  of  the  arrivals  as  possible,  two  widely 
varying  models  were  chosen:  water-plexiglas-plaster ,  and  water-glass- 
aluminum.  Primary  attention  will  be  directed  toward  the  water-plexiglas- 
plaster  model,  although  several  aspects  of  the  results  obtained  from 
the  water-glass-aluminum  model  will  be  discussed.  Because  of  the 
limiting  assumption  that  all  arrivals  experiencing  two  or  more  reflections 
from  the  solid-solid  interface  may  be  neglected,  it  should  not  be  expected 
that  all  experimental  arrivals  will  be  predicted  theoretically. 


(Text  Continued  on  Page  37) 
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FIGURE  5a.  EXPERIMENTAL  PRESSURE  RESPONSE  FOR  A  WATER- PLEXIGLAS- PLASTER  MODEL 
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FIGURE  5c.  EXPERIMENTAL  PRESSURE  RESPONSE  FOR  A  WATER-PLEXIGLAS-PLASTER  MODEL 
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FIGURE  5d.  EXPERIMENTAL  PRESSURE  RESPONSE  FOR  A  WATER-PLEXIGLAS-PLASTER  MODEL 
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FIGURE  5e,  EXPERIMENTAL  PRESSURE  RESPONSE  FOR  A  WATER-PLEXIGLAS-PLASTER  MODEL 
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FIGURE  6b.  EXPERIMENTAL  PRESSURE  RESPONSE  FOR  A  WATER-GLASS-ALUMINUM  MODEL 


WATER  -  PLEXIGTJVS  -  PUNSTER  MODEL 


Table  2  contains  the  computed  theoretical  travel  times  in  the  water- 
plexiglas-plaster  model  when  H  =  0.05  centimeters,  and  z  =  0.0  centi¬ 
meters.  These  theoretical  arrival  times  have  been  included  in  Figure  5 
so  that  they  might  easily  be  compared  with  the  experimental  results. 

It  has  been  shown  previously^  that  a  Rayleigh  wave  or  a  critically 
refracted  shear  wave  cannot  be  propagated  when  the  shear  velocity  of  the 
solid  at  the  fluid-solid  interface  is  less  than  the  compressional  sound 
velocity  of  the  fluid.  Hence,  since  such  a  relationship  between  elastic 
constants  is  present  in  the  water  -  plexiglas  -  plaster  model,  the 

and  Rayleigh  arrivals  do  not  exist.  It  can  be  shown  in  a  manner  similar 
to  that  utilized  above  that  the  ^1^2^2^^2^2^1 

arrivals  may  not  exist  in  this  model.  Also,  when  the  elastic  parameters 
of  the  two  solids  are  analyzed  with  respect  to  Scholte's  criteria  for  the 
existence  of  surface  waves  at  a  solid-solid  interface  (Figures  3  and  4), 
it  can  be  seen  that  such  an  arrival  is  not  possible  for  this  model. 

Referring  to  Figure  5a,  it  can  be  seen  that  for  this  range  of 
10  centimeters,  the  first  arrival  (38  microseconds)  appears  to  be  the 
wave,  the  next  arrival  (45  microseconds)  is  the  ^i^^2^3^2^1 

arrival,  and  the  third  arrival  (49  microseconds)  appears  to  be  the 
PlP2RP2^1  small  peaks  appear  at  58  and  72  microseconds 

which  have  not  been  identified.  At  69  microseconds  there  is  a  relatively 
strong  arrival  which  appears  to  be  a  superposition  of  the  P^P2RS2P|^> 

^l^^l  Following  the  direct  arrival 

there  is  a  region  of  small  oscillations  which  was  also  noted  in  the  fluid- 
homogeneous  solid  problem.  This  oscillation  is  again  related  to  the  fact 
that  we  have  an  explosive  type  point  source  instead  of  a  delta-excited 
line  source  as  the  theory  assumed.  The  direct  pulse  is  about  4  micro¬ 
seconds  long  and  the  reflected  (Pj^RP^)  pulse  occurs  soon  after  the 

beginning  of  the  direct  pulse.  The  Pj^RPj^  arrival  experiences  a  large 

phase  change  when  it  is  reflected  at  this  range,  so  strong  interference 
should  be  expected  in  this  region  of  the  pressure  response  curve  as  the 
arrival  times  for  the  direct  and  P^RP^  arrivals  are  almost  the  same. 

The  arrival  at  90  microseconds  is  due  to  the  P]^^2^3^2^1 

A  very  large  pressure  increase  begins  at  about  97  microseconds  and 
reaches  a  peak  at  about  100  microseconds.  Three  arrivals,  the  P^S^S^S^Pj^, 

^1^2*^^2^1*  Stoneley  (Liquid-Solid)  are  predicted  in  this  region. 
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TABLE  2 


THEORETICAL  ARRIVAL  TIMES  FOR  THE  WATER  -  PLEXIGLAS  -  PLASTER 
MODEL  (H  =  0.05  cm,  z  =  0.0  cm,  x  =  Variable) 


Wave  Type 

Arrival 

Time  in  Microseconds 

10.0 

X  = 

cm 

12.0 

X  = 

cm 

14.0 

X  = 

cm 

16.0 

X  = 

cm 

18.0 

X  = 

cm 

t(D) 

68.5 

82.2 

95.8 

109.5 

123.3 

t(PjRPp 

68.5 

82.2 

95.8 

109.5 

123.3 

t(PiPjPi) 

38.8 

46.4 

54.0 

61.6 

69.3 

t(P^S2P^) 

- 

- 

- 

- 

- 

46.8 

53.9 

60.9 

67.9 

75.2 

96.4 

107.5 

118.5 

129.5 

140.8 

88.7 

95.5 

102.5 

109.6 

116.8 

t(P^P2S3P2Pp 

- 

- 

" 

- 

- 

t(PjP2P3S2Pp 

67.6 

74.7 

81.7 

88.8 

96.0 

« 

" 

- 

- 

- 

t(PjP2RP2Pp 

48.9 

54.2 

61.2 

68.0 

76.0 

t(PjS2RS2Pp 

96.7 

109.4 

122.6 

135.9 

149.6 

t(P3P2RS2Pp 

67.5 

74.4 

81.2 

88.9 

96.3 

89.6 

97.0 

104.7 

112.4 

120.0 

- 

- 

- 

- 

- 

t (Rayleigh) 

- 

- 

- 

- 

- 

t (Stone ley :Liq, -Sol.) 

97.0 

116.2 

136.0 

155.5 

174.5 

t (Stone ley: Sol, -Sol.) 

- 

- 

- 

- 

- 
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Pressure  response  curves  obtained  at  larger  ranges  will  shov;  that  this 
large  peak  is  primarily  due  to  a  Stoneley  (Liquid-Solid)  arrival* 

Many  of  the  arrivals  occurring  after  the  Stoneley  peak  are  due  to  re¬ 
flections  disallowed  in  the  theoretical  derivation.  It  will  be  shown 
that  the  large  arrival  at  about  142  microseconds  is  also  a  Stoneley 
(Liquid- Solid)  arrival. 

The  first  three  arrivals  in  Figure  5b  (range  =  12  centimeters) 
again  appear  to  be  the  P^P^P^,  ^1^2^2^  arrivals, 

occurring  at  47,  52,  and  56  microseconds  respectively.  The  two  small 
arrivals  at  63  and  68  microseconds  have  not  been  identified.  The  peak 
at  75  microseconds  is  a  superposition  of  the  P^P^RS^P^  and  ^]^^2^3^2“1 

arrivals.  The  direct  and  reflected  (P^RP^)  arrivals  occur  at  about 

82  microseconds,  followed  by  the  usual  region  of  oscillations,  seme  of 
which  are  relatively  large  in  this  pressure  response  curve.  The 
^1^2^2^^2^2^1  ^1^2^3^2^1  predicted  at  97  microseconds, 

but  appear  to  be  somewhat  distorted.  The  arrival  at  102  microseconds 
is  unexplained.  The  large  arrival  with  a  peak  at  108  microseconds  is 
predicted  to  be  the  ^]^^2^3^2^1  ^  weak  arrival,  possibly  the 

P^S^RS^Pj^  occurs  at  112  microseconds.  It  must  be  noted  that  many  of 

the  reflected  arrivals  have  probably  experienced  phase  changes  greater 
than  90  degrees;  and  hence,  these  arrivals  may  be  negative  rather  than 
positive  peaks.  The  large  Stoneley  (Liquid-Solid)  arrival  begins  at 
114  microseconds  and  reaches  a  peak  at  117  microseconds.  The  strong 
arrival  at  about  156  microseconds  will  be  identified  as  a  Stoneley 
(Liquid-Solid)  arrival  later. 

Referring  to  Figure  5c,  it  is  now  seen  that  the  two  small  arrivals 
at  63  and  68  microseconds  are  unidentified.  Also,  two  moderate  arrivals 
at  127  and  130  microseconds  are  unidentified.  Of  course,  the  two  pre¬ 
dominant  waveforms  are  the  Stoneley  arrival  at  135  microseconds  and 
the  other  large  arrival  at  175  microseconds  which  will  be  discussed 
shortly. 

Figures  5d  and  5e  show  that  the  Stoneley  arrival  continues  to  be 
the  predominant  waveform,  although  the  ^j^^2^3^2^1  beginning 

to  show  a  rather  large  amplitude.  Several  arrivals  continue  to  be 
unidentified. 


Figure  7  shows  the  travel  time  versus  range  curves  for  the  Stoneley 
(Liquid-Solid)  arrival  and  for  the  large  unidentified  arrival  that 
usually  occurs  about  40  microseconds  after  the  Stoneley  arrival.  The 
inverse  slopes  of  these  curves  represent  the  velocity  with  which  the 
waves  are  traveling.  The  theoretical  Stoneley  velocity  is  1.03  x  10^ 
centimeters  per  seconds,  while  the  experimental  velocity  is 
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FIGURE  7.  TRAVEL  TIME  VERSUS  RANGE  CURVES  FOR  THE  STONELEY  ARRIVALS 


1.09  X  10  centimeters  per  second.  The  large  unidentified  arrival  that 
is  suspected  to  be  a  Stoneley  arrival  has  a  velocity  of  1.11  x  10^ 
centimeters  per  second.  Since  this  velocity  is  considerably  less  than 
the  sound  velocity  in  either  of  the  solids,  it  is  strongly  indicated 
that  this  arrival  must  be  identified  as  a  Stoneley  arrival.  But,  the 
Stoneley  wave  is  propagated  only  along  the  interface  (i.e.,  horizontally) 
so  this  latter  arrival  cannot  represent  a  reflection  of  the  Stoneley  wave. 
The  possibility  remains  that  a  compressional  wave  in  the  fluid  could  be 
refracted  as  it  enters  the  upper  solid,  reflected  at  the  solid-solid 
interface,  and  then  strike  the  fluid-solid  interface  at  such  an  angle  so 
that  it  would  force  the  formation  of  a  Stoneley- type  arrival.  But,  if 
this  is  the  situation,  geometrical  considerations  show  that  such  an  ar¬ 
rival  could  not  occur  in  this  model  unless  the  range  was  greater  than 
15  centimeters.  Thus,  the  exact  reason  for  the  existence  of  the  latter 
"Stoneley- like"  arrival  is  unknown.  It  is  apparent  that  further  research 
related  to  the  rather  strange  physical  mechanisms  by  which  surface  waves 
are  formed  and  propagated  will  be  necessary  to  solve  this  problem. 


WATER  -  GLASS  -  ALUMINUM  MODEL 

The  theoretical  travel  times  in  the  water  -  glass  -  aluminum  model 
when  H  =  0.5  centimeters  and  z  =  0,0  centimeters  are  given  in  Table  3. 

As  before,  these  theoretical  times  have  been  included  in  Figure  6  so 
that  they  might  easily  be  compared  with  the  experimental  results. 

The  first  arrival  in  Figure  6a  appears  to  be  due  to  both  the 
P^P2P^and  Pj^P2P2P2^1  next  two  arrivals  are  due  to  the 

^1^2^ 3^2^ I  about  35  microseconds,  the  direct 

arrival  occurs,  supported  by  the  Rayleigh,  and  ^^^52?^ 

arrivals.  The  reflected  (P^RP^)  arrival  does  not  now  interfere  with  the 

direct  arrival  as  much,  as  there  is  almost  a  1  microsecond  difference  in 
arrival  times.  The  peak  occurring  at  about  41  microseconds  is  the 
Stoneley  (Liquid-Solid)  arrival.  No  Stoneley  (Solid- Solid)  arrival 
can  exist  in  this  model.  In  addition,  several  of  the  reflection  and 
refraction  arrivals  cannot  occur,  as  can  be  seen  in  Table  3. 

Referring  to  Figure  6b,  the  first  arrival  (26  microseconds)  is  now 
the  Th^i  is,  the  arrival  that  is  critically  re¬ 

fracted  at  the  solid-solid  interface  has  a  shorter  travel  time  than 
the  arrival  that  is  critically  refracted  at  the  fluid-solid  interface. 
The  arrival  at  30  microseconds  is  evidently  due  to  a  superposition  of 
the  P^P2Fj^  and  P^P2RP2P^  arrivals.  The  arrival  is  predicted 

to  occur  31.6  microseconds,  but  no  such  arrival  is  present.  The  con¬ 
tribution  at  about  t  =  36  microseconds  is  due  to  the  P^P2RS2P^and 
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TABLE  3 


THEORETICAL  ARRIVAL  TIMES  FOR  THE  WATER  -  GLASS  -  ALUMINUM 
MODEL  (H  =  0.5  cm,  z  =  0.0  cm,  x  =  Variable) 


Wave  Type 

Arrival  Time 
X  =  5,0  cm 

in  Microseconds 

X  =  7 . 5  cm 

t(D) 

34.2 

51.3 

t(P^RP^) 

35.0 

52.3 

21.3 

28.  9 

tCPjS^Pj^) 

34.4 

50,4 

t(P^P2P3P2Pp 

22.1 

26.9 

t(P^S2S3S2Pp 

- 

- 

31.6 

36.4 

t(P^P2S3P2P3) 

- 

- 

26.8 

31.6 

t(P^P2S3S2P^) 

- 

- 

t(P^P2RP2Pp 

23.3 

30.3 

t(Pj^S2RS2P3) 

40.0 

54.0 

t(P^P2RS2P3) 

29.1 

36.5 

t(P^P2S2RS2P2Pj^) 

- 

- 

t(P^S2P2RP2S2Pj^) 

35.7 

43.2 

t (Rayleigh) 

35.0 

52.4 

t (Stone ley :Liq. -Sol. ) 

40.0 

60.0 

t(Stoneley:Sol. -Sol.) 

- 

- 
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*  The  P^S2p2R^2^2^1  occurs  at  42,5  micro¬ 

seconds,  although  it  is  quite  weak.  Also,  a  small  disturbance  at  about 
40  microseconds  is  believed  to  be  the  arrival.  Again,  several 

arrivals  support  the  direct  pulse  at  about  51  microseconds.  The  Stoneley 
zero  crossing  is  at  58,8  microseconds,  while  the  Stoneley  peak  occurs 
at  about  60  microseconds. 


CONCLUSIONS 


The  most  important  conclusions  that  can  be  stated  at  this  time  are: 

1.  The  theoretical  approach  utilized  appears  to  be  quite  ap¬ 
plicable  to  problems  of  the  type  considered  here.  Arrival  time  pre¬ 
dictions  can  be  obtained  with  little  effort,  and  exact  amplitude  data 
should  be  available  when  the  mathematical  development  can  be  completed. 

2.  The  extended  theory  developed  in  this  report  is  accurate  in 
predicting  most  expected  arrivals  at  short  ranges  in  the  experimental 
model,  although  several  experimental  arrivals  are  received  that  are  not 
predicted.  Further  development  of  the  theory,  which  will  give  waveform 
data,  should  prove  useful  in  determining  the  origin  of  these  unexpected 
arrivals , 

3.  Except  for  rare  situations,  Stoneley  waves  cannot  exist  at 
real  solid-solid  interfaces.  It  must  be  emphasized  that  such  waves  are 
theoretically  possible  under  certain  conditions,  but  that  these  con¬ 
ditions  are  so  stringent  that  they  can  be  rarely  satisfied. 

4.  An  additLonal  Stoneley- type  wave  appears  to  exist  at  the  fluid- 
solid  interface  because  of  a  reflection  from  the  solid-solid  interface, 
although  this  cannot  be  substantiated  by  a  geometrical  analysis, 

5.  The  work  performed  on  this  phase  of  the  project  continues  to 
illustrate  the  large  errors  in  wave  propagation  predictions  near  an 
interface  that  can  occur  when  a  hard  bottom  is  theoretically  considered 
to  be  a  fluid, 

6.  Although  this  paper  has  considered  only  a  fluid-two  layer- 
solid  model,  the  results  are  useful  when  considering  other  models.  The 
principal  contribution  to  the  pressure  response  waveform  obtained  in 
the  fluid  near  a  hard  bottom  will  often  be  the  Stoneley  wave  propagated 
at  the  fluid-solid  interface,  although  several  such  waves  may  be  present 
if  the  bottom  contains  a  number  of  layers. 
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DIRECTION  OF  FUTURE  WORK 


It  is  apparent  that  much  research  remains  to  be  done  in  the  field 
of  elastic  wave  propagation  near  an  interface.  Many  questions  remain 
unanswered,  as  has  been  noted  throughout  this  report.  Certainly,  it 
seems  that  surface  waves  must  be  responsible  for  a  major  part  of  the 
pressure  response  when  one  considers  the  propagation  of  low-frequency 
sound  near  a  hard  ocean  floor,  and  hence  a  thorough  understanding  of 
such  wave  propagation  is  important. 

The  problem  of  elastic  wave  propagation  has  been  under  investi¬ 
gation  since  before  the  turn  of  the  century,  and  research  in  this  field 
will  certainly  continue  well  into  the  future.  Because  of  the  complexity 
of  the  mathematics  involved  in  the  rigorous  solution  of  any  elastic 
wave  propagation  problem,  many  limiting  assumptions  are  almost  always 
necessary.  Often  it  is  necessary  to  accept  assumptions  that  make  ex¬ 
perimental  verification  of  theoretical  results  difficult,  if  not  im¬ 
possible,  For  this  reason,  the  efforts  expended  in  this  field  at  this 
Laboratory  have  been  planned  so  that  experimental  and  theoretical  de¬ 
velopments  go  side-by-side.  Models  are  chosen  so  the  theoretical 
assumptions  necessary  to  solve  a  specific  problem  may  be  experimentally 
duplicated. 

Future  endeavors  on  this  project  should  include: 

1,  A  completion  of  the  mathematical  development  presented  in  this 
report--This  will  consist  of  performing  the  transformation  from  the 
variable  q  to  real  time,  and  then  programming  the  results  through  use 
of  a  computer  so  that  pressure  response  versus  time  curves  can  be 
easily  obtained  for  any  system  to  be  considered. 

2,  Inclusion  of  the  contributions  due  to  reflections  from  the 
water-air  interface--This  will  be  a  mathematical  and  experimental 
development  quite  similar  to  that  presented  in  this  report. 

3,  Comparison  of  theoretical  and  model  work  with  data  obtained 
over  a  soft  bottom  for  an  explosive  type  sound  source--This  data  is  now 
available  and  preliminary  analysis  should  begin  in  the  near  future. 
Although  most  of  the  model  work  has  been  done  with  rigid  materials,  it 
must  be  stressed  that  the  theory  is  also  applicable  to  softer  bottom 
materials.  Figure  8a,  b  and  c  contains  pressure  response  curves  for 
systems  consisting  of  (a)  recent  sediments,  (b)  consolidated  sediments, 
and  (c)  metamorphic  rock.  It  can  be  seen  that  Stoneley  wave  contri¬ 
butions  are  substantial  even  for  the  recent  sediment  bottom.  Since  the 
Stoneley  wave  is  symmetrical,  most  of  the  response  following  the  direct 
arrival  can  be  attributed  to  it.  Hence,  an  analysis  of  the  available 
data  should  show  surface  wave  contributions. 
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FIGURE  8a.  THEORETICAL  PRESSURE  RESPONSE  CURVE  FOR  A 
WATER  -  RECENT  SEDIMENT  SYSTEM 


FIGURE  8b.  THEORETICAL  PRESSURE  RESPONSE  CURVE  FOR  A  WATER 
CONSOLIDATED  SEDIMENT  SYSTEM 
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4.  Extension  of  the  theory  to  repetitive  impulsive  8ource8--Thi8 
would  consist  of  both  theoretical  and  model  work. 

5.  Consideration  of  a  harmonic  type  source — If  possible,  a  **point** 
harmonic  type  sound  source  should  be  obtained  and  model  work  done  using 
this  source.  The  theory  considered  to  data  should  be  useful  with  suitable 
modification,  and  an  experimental  endeavor  will  prove  useful. 

6.  Investigation  of  the  scale  laws  governing  the  various  types 
of  waves  present  in  the  model  work--Thi8  would  include  determination  of 
amplitude  decay  rates  and  frequency  content  changes  with  varying  dimensional 
parameters. 
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APPENDIX  A 


DERIVATION  OF  THE  REFLECTION  COEFFICIENTS  A(q)  AND  N(q) 


Equations  (1)  and  (2)  define  the  11  displacenient  potentials  con¬ 
sidered  in  this  problem.  These  potentials  are  subject  to  the  following 
boundary  conditions  (See  Figure  2)  : 

1.  The  normal  stresses  (azz)  are  continuous  at  z  =  H  and  z  =  H  +  J 

2.  The  normal  displacements  (w)  are  continuous  z  =  H  and  z  =  H  +  J 

3.  The  tangential  stresses  (axz)  are  zero  at  z  =  H  and  continuous 

at  z  =  H  +  J, 

4.  The  tangential  displacements  (u)  are  continuous  at  z  =  H  +  J. 

These  boundary  conditions  give  10  boundary  equations  which  can  be 
used  to  solve  for  the  10  quantities  A,  B,  C,  D,  E,  F,  G,  I,  M  and  N. 
Hence,  instead  of  the  3  equations  and  3  unknowns  in  the  fluid-single 
solid  problem,  we  now  have  10  equations  and  10  unknowns.  If  we  had  not 
limited  the  problem  to  a  maximum  of  two  reflections  from  the  solid-solid 
interface,  even  more  potentials,  and  hence  unknowns,  would  have  been 
present . 

Since  the  present  problem  concerns  the  response  in  the  fluid,  it  is 
only  necessary  to  determine  A(q)  and  N(q). 

The  displacements  and  stresses  are  given  by^: 


zz 


f-  ^2s  ^2,s^ 

.  pf  (1  -  2  ij)  ^  +  2b"  (4-f-  - 

^  a  dt  dz  3x9z 


r:,2  s  „  2  s  ,  s 

p  '  -  +  2b" 


r  ir  :»2  f 

f  cc  X 

a  =  P  ("2  ^ 
ot 


zz 


a  "  *  0 

xz 
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s 


w 


^  _  Ml  ^ 

bz  dz  '  ^  bz 


s  ^  .  Ml 

Applying  the  boundary  conditions  and  inserting  the  displacement 
potentials  (2a  -  j)  into  the  boundary  equations,  we  obtain: 

-  2b V)]  ^  •{ap^lb^qb-.'^'’''*] 

-  -p‘e-“™ 

+  B[p'e-“»’“]  -  c[,b-“’«]  .  rr‘“™  (A-2) 

B[2b^.'qe"‘''“*''‘]  +  c[(l  -  2b^p^)e'‘“'’’“]  -  0  (A-3) 

b[p  ^(1  -  2b^q^)e'‘“*'*"  -  c[2p*  Vb 'qe'*'”*’’ 

+  Il[p‘l(l  -  2bV)p“'““'<"  *  ■'’]  +  E[2p*Vb'qe-‘»'’'<»  +  ■’>] 

+  l[2p  Vqk'e"‘'“'‘'*‘'  -  h[p  ^(1  -  2k^q^)e'“*’ -  0  (A-A) 

l(2p‘lb^'qp-‘'"'<"  *  c[p“l(l  -  2bV)b'‘'“'’''''  ^  ■’>] 

-  ll[2p‘V.’q.-‘“'<“  +  ■'>]  A  e[p'\i  -  2bV)b-“'‘'*  *  ■’>] 

-  l[p  ^(1  -  2k^q^)e'‘“'‘'*"  *  '^5  -  'i}”  ^E^E'qp'^”®  *'*  0  (A-5) 
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-  +  J>]  -  c[b'e-“'>'<“  +  J)]  -  d[, 

+  E[b'e-‘»'=’<»  ■»]  +  +  Jj]  ^  +  J)  j  .  „ 

-  B[a'e-‘”‘‘'<''  +  J)]  +  +  ''>]  +  DTa’a-l”*' <»  +  ■'>1 

+  +  ■'>]  +  +  J)]  .  i[qa-“l>’("  +  J)]  .  0  <A-7) 

d[p°1(1  -  2b%^)a-“»'(«  +  2-’)]  +  E[2p'V,b'a-*’<“  ^ 

+  ffp'Nl  -  2bV)q-''“'<"  ^  '">]  -  G[2p"V,b'a-“'»  +  2J)] 

.^pf^-l».r(K  +  2J)T  .  „ 

e[(1  -  2bV)q'“’’'<'‘  "■'']  +  F[2b2a',a-“‘'<“  +  “>] 

+  g[(1  -  2b\^)e'‘“'’'<"  ■*■  -  D[2b^«'qe'““'^'*  ■*■  -  0  (A-9) 

-  +  2J)]  .  E[,,-Wb’(H  +  2J)J  ^  4-  2J)J 

-  ((qa-^X''"  +  4  .(ra-^“^<“  *  -  0  (A- 10) 

This  set  of  10  equations  can  be  written  very  compactly  in  matrix 

form: 
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It  is  now  theoretically  possible  to  solve  for  each  of  the 
10  unknown  quantities  by  evaluating  two  10  by  10  determinants.  It  should 
be  noted  that  these  determinants  immediately  reduce  to  7  by  7  determinants 
due  to  the  fact  that  A(q) ,  B(q),  and  C(q)  are  not  dependent  on  the  other 
unknown  quantities.  Such  evaluations  will  yield  equations  (3)  and  (4). 


(Reverse  Page  34  Blank) 
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APPENDIX  B 


THE  REFRACTION  AND  REFLECTION  ARRIVALS 


It  has  been  shown^*^  that  in  an  infinite  fluid  only  a  compressional 
type  of  wave  motion  can  occur,  while  in  an  infinite  solid  both  com¬ 
pressional  and  shear  waves  are  possible.  Considering  the  liquid-solid 
interface  illustrated  in  Figure  Bl,  it  can  be  shown  that  incident  com¬ 
pressional  waves  from  the  liquid  can  produce  compressional  and  shear  waves 
in  the  solid,  and  also  that  a  surface  wave  which  is  tied  to  the  interface 
can  occur.  For  certain  angles  of  incidence  for  the  impinging  energy, 
critically  refracted  waves  can  occur.  Furthermore,  reflections  from  the 
interface  are  possible.  Hence,  for  this  fluid-solid  problem,  one  re¬ 
flection  arrival,  two  refraction  arrivals,  one  direct  arrival,  and  surface 
arrivals  are  predicted. 


P 


FIGURE  Bl.  THE  LIQUID  -  SOLID  INTERFACE 


The  problem  under  consideration  in  this  report  concerns  a  fluid  overlying 
a  two-layer  solid  bottom.  The  prediction  of  the  various  reflections  and 
refraction  arrivals  now  becomes  much  more  difficult.  In  fact,  as  dis¬ 
cussed  previously,  there  are  an  infinite  number  of  wave  arrivals  unless  a 
limiting  condition  is  applied  to  the  problem.  In  the  mathematical 
derivation  section  of  this  report  the  assumption  was  made  that  any  waves 
existing  after  two  reflections  from  the  lower  (solid-solid)  interface 
could  be  neglected.  As  stated  previously,  this  is  not  always  a  valid 
assumption,  but  for  the  purpose  of  this  investigation  it  is  acceptable. 
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Figure  B2,  B3,  B4,  and  B5  show  the  ten  refraction  arrivals  that  must 
be  considered  in  the  problem. 


FIGURE  B2,  THE  ^1^2^3^2^r  ^1^2^3^2^1 

REFRACTION  ARRIVALS 


If -  X  - H 


FIGURE  B3.  THE  ^1^2^3^2^1  ^^^^^^TION  ARRIVALS 


(Text  Continued  on  Page  58) 
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FIGURE  B5.  THE  P^S2P3P2Pj^  AND  P^S2S^P2P^  REFRACTION  ARRIVAL 
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Figures  B6  and  B7  illustrate  the  five  reflection  arrivals  that  must 
be  considered: 


|| -  X  - ^ 

Source 


FIGURE  B6.  THE  RE^'LECTION  ARRIVALS 


|« -  X - ^ 

Source 


FIGURE  B7.  THE  P^P2RS2P^  AND  P^S2RP2P^  REFLECTION  ARRIVALS 
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P  .8,k 


FIGURE  B9.  THE  P^S2P2RP2S2^1  ^^^FRACTION  -  REFLECTION  ARRIVAL 


Also, we  should  expect  surface  or  interface  waves  at  both  the  liquid- 
solid  and  solid-solid  interfaces,  as  has  been  shown  previously.  Hence, 
for  the  assumed  liquid- two  solid  problem,  we  should  expect  one  direct 
arrival,  ten  refraction  arrivals,  five  reflection  arrivals,  two  refraction 
reflection  arrivals,  and  the  surface  arrivals  from  both  interfaces.  In 
other  words,  there  are  now  over  twenty  arrivals  to  analyze,  rather  than 
the  six  arrivals  found  in  the  liquid-solid  problem. 


It  has  been  shown  previously^  that  the  arrival  times  for  the 
and  refraction  arrivals  are  given  by: 


t 


P  P  P 
1  2^1 


and 


t 


(B-l) 


(B-2) 


The  derivations  of  the  arrival  tiroes  of  the  remaining  critically 
refracted  arrivals  are  similar,  so  only  the  will  be 

derived  and  it  will  be  used  as  a  guide  to  obtain  the  other  arrival  tiroes. 
The  symbols  used  are  defined  in  Figure  BIO. 


FIGURE  BIO.  DEFINITION  OF  SYMBOLS  FOR  THE  CALCULATION  OF  THE 
OF  THE 
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From  Snell's  law, 


,  Q  sin  0  sin  9 

*^°  °1 _ ^ 

2  a  ’  a  “  g 

From  Figure  BIO: 

(1^)^  =  (H)^  + 

(§2)^  =  (H  -  z)^  +  (x^)^ 


(53)^  =  (J)^  +  (X3)^ 


X  x^ 
sin  e  =  ^  =  — 
^1  ^2 


sin  0  =  T~ 

C2  ?3 

Therefore,  the  travel  times  in  the  fluid  are  given  by 

.  ^  il  ^  fiH  . 

1  c  2  2^ 

c(g  -  c^) 

and 


,  _  12  ^  g(H  -  2) 

4  -  c  ,  2  4  • 

c(g  -  C  ) 

Likewise,  the  travel  times  in  the  upper  solid  are  given  by 


Jg 


*^3  “  a  ~  ,  j  \ 

a(g  -  a^) 
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The  travel  time  in  the  lower  solid  is  given  by 


X  -  -  2X3  ^  -  Z)c  .  2aJ 

A  g  g  2  2:  ^,2  2.^ 

g(g  -  c  )  g(g  -  a  ) 

Therefore,  the  total  travel  time  for  the  arrival 

is  given  by 


P  P  P  P  P 
1  2^3  21 


2  +  (2H  -  z) 
8  eg 


(g 


(B-3) 


The  travel  time  expressions  for  the  other  critically  refracted  arrivals 


can  be  found  in  a  manner  similar  to  that  used  to  find  t 


travel  times  are  given  by: 


X  .  (2H  -  z)  ..  2  2^  .  to 

t„  „  „  „  „  =  -  +  (k  -  c  )  +  —  (k  -  b  ) 


P  P  P  P  P 
1  2^3  2  1 


2J  ,.  2  ,  2, 

bk 


^"2Wl 


X  .  r2H  -  ,  2  2.^  .  2J  .  2  .2. 

=  7  +  — 7:;^ —  (g  ■  )  +  —  (g  ■  ) 

g  eg  bg 

*  7  +  +  S  (k^  - 


ak 


,  _  X  (2H  -  z)(g^  -  c^ 

*’l^2^3^2^1  ^1^2^3^2^1  8  gc 


These 

(A-A) 
(A- 5) 
(A-6) 


+  -^  (g^  -  +  £  <8^  - 

ag  bg 


,  ^  X  ■  I2H  -  (k^  -  c^) 

^1^2SVl  ^^2SVl 


T  2  2^  T2  2^ 


(A-7) 


(A-8) 


The  travel  time  expressions  for  the  reflection  arrivals  are  not  as 
simple  as  the  equations  obtained  for  the  refraction  arrivals.  The  geo¬ 
metrical  approach  necessary  to  determine  the  travel  times  of  the  reflection 
arrivals  consists  of  an  iterative  process.  For  a  given  x,  z,  H,  J,  c,  a, 
and  b,  a  travel  path  is  found  that  satisfies  Snell’s  law  at  each  boundary 
and  still  passes  through  the  source  and  detector.  Then,  once  this  path  is 
known,  it  is  possible  to  compute  the  travel  time. 
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It  is  possible  to  obtain  an  algebraic  expression  for  the  travel 
times  of  the  refraction  reflection  arrivals.  Figure  Bll  defines  the 
symbols  that  will  be  used  in  the  travel  time  derivation  for  the 


FIGURE  Bll.  DEFINITION  OF  SYMBOLS  FOR  THE  CALCULATION  OF  THE  ARRIVAL 
TIME  OF  THE  ? ^RS ^  ARRIVAL 
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Snell's  law  states  that: 


sin  9 


sin  0. 


1 

a 


From  Figure  Bll: 


=  xj  + 


,2  2  ^  2 


2  2  2 
=  X5  +  (H  -  z)^ 


Xi  X 

sin  9  =  —  =  ■^ 

*=1  ^4 


sin  0  =  ^ 

3 

Therefore,  the  travel  times  in  the  fluid  are  given  by 


and 


*^1  c 


aH 


^  2  2  ' 
c(a  -  c  ) 


_  12  _  a(H  -  z) 

5  ■ '  ■  2  2  i 

c(a  -  c  ) 


The  travel  times  in  the  solid  are  given  by 

3  b  o  0  ^ 

b(a  -  b'^) 


and 


(*2  *4^  X  c(2H  -  z)  2bJ 

*^2  ^ - 1 -  “  -  “  ^  ■ 


®  /  2  2.^  .2  .2.^ 
a(a  -  c  )  a(a  -  b  ) 
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Thus,  the  total  travel  time  for  the  given  by 


t 


X  .  (2H  -  z)  .  2  2.  ^  2J  .2 

^  - -  (a  -  c  )  +  "-r  (a 

a  ac  ab 

Likewise,  it  can  be  shown  that 


b^) 


^P^s^PzRPaSzPi 


X  .  (2H  -  2j 
b  be 


(b' 


(b 


(A-9 


(A- 10) 
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APPENDIX  C 


EXPERIMENTAL  INSTRUMENTATION 


The  experimental  apparatus  employed  in  this  project  was  similar  to 
that  used  by  Roever  and  Vining  in  their  original  verification  of  Strick*s 
theory,^  A  detailed  description  of  the  apparatus  is  given  in  Appendix  D 
of  Reference  1.  Only  a  brief  description  of  the  apparatus  will  be  in¬ 
cluded  here,  along  with  details  of  the  few  modifications  necessary  in 
the  apparatus  to  investigate  the  stratified  bottom  problem. 

The  spark  source  and  detection  system  remains  unchanged  from  the 
previous  work.  The  spark  source  consisted  of  a  coaxial  spark  gap  to 
which  was  applied  the  power  from  a  1  |xf  capacitor  charged  to  812  volts. 
The  peak  current  during  the  discharge  was  over  1000  amperes,  and  the  dis¬ 
charge  lasted  for  about  4  microseconds.  The  size  of  the  spark  gap  was 
about  1  millimeter,  so  it  could  be  considered  to  be  a  point  source. 


The  detector  was  a  Glennite  UP-800c  ultrasonic  probe.  This  is 
simply  a  commercial  barium  titanate  detector  designed  for  minimal  dis¬ 
turbance  of  the  sound  field.  The  sensitivity  was  -153(tl)db  (ref: 
Iv/dyne/cm^)  between  10  kc  and  400  kc,  with  good  response  down  to  2  kc 
and  up  to  2  me. 


The  probe  output  led  to  a  Tektronic  545  oscilloscope  which  was 
triggered  by  a  low-resistance  shunt  in  the  discharge  circuit.  Time 
resolution  was  ±1  microsecond,  and  was  limited  by  the  finite  size  of  the 
source  and  detector. 


A  plaster  block  (12**  x  12*'  x  10'*)  used  in  the  previous  work  was 
also  employed  in  this  problem  as  one  of  the  materials  used  for  the 
lower  semi-infinite  solid.  A  large,  solid  aluminum  cylinder  was  also 
used  for  this  purpose.  The  layer  materials  used  were  a  1.27  cm  slab 
of  glass  and  a  3.81  cm  slab  of  plexiglas.  The  elastic  parameters  of 
these  materials  are  given  in  Table  1, 

The  question  of  bonding  at  the  solid  interfaces  had  to  be  considered 
carefully.  In  theoretical  derivations  the  solid-solid  interface  is  con¬ 
sidered  to  be  **welded.**  That  is,  it  is  assumed  that  the  tangential 
stresses  and  displacements  are  continuous.  Hence,  should  the  solid 
layers  be  bonded  together,  or  is  it  sufficient  to  simply  place  them 
together?  A  search  of  the  available  literature  showed  about  an  equal 
division  between  the  two  techniques,  A  simple  solution  to  the  bonding 
difficulty  was  as  follows: 

First,  data  was  obtained  for  the  water  -  glass-  aluminum  model  with 
the  glass  simply  lying  on  the  aluminum  block,  x  was  6.5  cm,  H  was  0.2  cm. 
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and  z  was  0,0  cm.  Then,  the  glass  and  aluminum  were  bonded  with  Barge 
cement  and  the  experiment  was  repeated.  The  two  sets  of  data  are  com* 
pared  in  Figure  Cl.  It  can  be  seen  that  very  little,  if  any,  change 
occurred  because  of  the  bonding  of  the  two  materials.  Hence,  it  was 
concluded  that  bonding  was  not  necessary  in  this  problem,  although  if 
the  layer  surfaces  were  not  smooth  and  flat  the  bonding  would  probably 
have  been  of  considerable  importance. 
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FIGURE  Cl,  ILLUSTRATION  OF  THE  EFFECT  OF  BONDING  THE 
GLASS  AND  ALUMINUM  LAYERS 
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foiBi  of  tKo  oxtondod  thoory  ia  not  obtoinod  in  cloaod  form^  it  com  pradict  oil  form  of  tha  axtandad  thaory  is  not  obtoinad  in  closad  form,  it  com  pradict  all 

vrlirala  at  abort  rangaa.  It  ia  concludad  that  tha  axiatanca  of  Stonalay  wavaa  arrlvaia  at  short  ranqaa.  It  la  concludad  that  tha  aziatanca  of  Stonalay  wavas 

at  a  aoiicbaoiid  intarfoca  ia  uniihalyf  although  raflactiona  from  tha  aolid-aoUd  at  a  aoUd-solld  intarfaca  is  unlikaly,  although  raflactions  from  tha  soUd^aoUd 

intarfaca  moy  ba  raaponaibla  for  tha  fonnation  of  additional  Stonalay  wavaa  ot  intarfaca  may  ba  rasponsibla  for  tha  formation  of  odditionoU  Stonalay  wavaa  at 

tha  flui^aoUd  intarfaca.  tha  fluid- solid  intarfaca. 
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